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Abstract 

The geometrical diffraction theory, in the sense of Keller, is here reconsidered as an obstacle 
problem in the Riemannian geometry. The first result is the proof of the existence and the 
analysis of the main properties of the diffracted rays, which follow from the non-uniqueness of 
the Cauchy problem for geodesies in a Riemannian manifold with boundary. Then, the axial 
caustic is here regarded as a conjugate locus, in the sense of the Riemannian geometry, and the 
results of the Morse theory can be applied. The methods of the algebraic topology allow us 
to introduce the homotopy classes of diffracted rays. These geometrical results are related to 
the asymptotic approximations of a solution of a boundary value problem for the reduced wave 
equation. In particular, we connect the results of the Morse theory to the Maslov construction, 
which is used to obtain the uniformization of the asymptotic approximations. Then, the border 
of the diffracting body is the envelope of the diffracted rays and, instead of the standard saddle 
point method, use is made of the procedure of Chester, Friedman and Ursell to derive the 
damping factors associated with the rays which propagate along the boundary. Finally, the 
amplitude of the diffracted rays when the diffracting body is an opaque sphere is explicitly 
calculated. 
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1 Introduction 

Classical geometrical optics fails to explain the phenomenon of diffraction: the existence of non- 
zero fields in the geometrical shadow. In several papers |l], ||, [| Keller proposed an extension of 
classical geometrical optics to include diffraction (see also the book by Bouche, Molinet and Mittra 
and references therein, where all these results have been collected and clearly exposed). This 
modification basically consists in introducing new rays, called diffracted rays, which account for 
the appearance of the light in the shadow. The most clear and classical example of diffracted rays 
production is when a ray grazes a boundary surface: the ray splits in two, one part keeps going as 
an ordinary ray, whereas the other part travels along the surface. At every point along its path this 
ray splits in two again: one part proceeds along the surface, and the other one leaves the surface 
along the tangent to the surface itself (see Fig. [l]) . Keller gives also a heuristic proof of the existence 
of these diffracted rays which is based on an extension of the Fermat's principle Q. In spite of 
these efforts the concept of diffracted rays still remains partially based on physical intuition. The 
first aim of this paper is to put on firm geometrical grounds the existence and the properties of the 
diffracted rays when the diffracting body is a smooth, convex and opaque object. To this purpose the 
diffraction problem is here reconsidered as a Riemannian obstacle problem, and then the diffracted 
rays arise as a consequence of the non-uniqueness of the Cauchy problem for the geodesies at the 
boundary of the obstacle |5], |(| (i.e. the diffracting body). Next, we are faced with the problem 
related to the caustic 'M, which is composed of the obstacle boundary and of its axis (axial caustic). 
In particular, since the latter can be regarded as a conjugate locus in the sense of the differential 
geometry, all the classical results of the Morse theory j|] must be formulated by taking into account 
the main geometrical peculiarity of the problem: the manifold we are considering has boundary. 
Finally, the rays bending around the obstacle can be separated in various homotopy classes by the 
use of the classical tools of algebraic topology. All these geometrical questions will be analyzed in 
Section ||. 

Classical geometrical optics corresponds to the leading term of an asymptotic expansion of a 
solution of a boundary value problem for the reduced wave equation. This term, which is generally 
derived by the use of the stationary phase method, gives approximations which have only a local 
character: they are not uniform. In particular, these approximations fail on the caustic; then 
a problem arises: how to patch up local solutions across the axial caustic. It is well known in 
optics that after a ray crosses a caustic there is a phase shift of —tt/2. It will be shown that this 
result can be derived in a very natural way by the use of the Maslov construction || 1C], which 



effectively allows for linking the patchwork occurring when the ray crosses the axial caustic and the 
geometrical analysis of Section |2[ 

The surface of the diffracting body is the envelope of the diffracted rays: it is a caustic and, 
consequently, the classical method of the stationary phase fails on it. By using a modified method 



of the stationary phase, due to Chester, Friedman and Ursell [11], we derive a countably infinite 
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Figure 1: Geometric representation of the propagation of the creeping waves into the shadow of a 
spherical obstacle. 

set of factors which describe the damping of the creeping waves along the surface: these damping 
factors depend on the obstacle curvature. These results, which lead to the Ludwig-Kravtsov |L2], [13|] 
uniform expansion at the caustic, are briefly discussed in subsection |3.2| . 
From these introductory considerations the following geometrical ingredients emerge: 

i) the non-uniqueness of the Cauchy problem at the boundary of the manifold in a Riemannian 
obstacle problem; 

ii) the correspondence between the homotopy classes of the fundamental group of the circle and 
the number of crossings through the axial caustic; 

hi) the Maslov phase-shift associated with the crossing-number of the ray; 

iv) the relationship between the curvature of the obstacle and the damping of the creeping waves 
along the surface of the obstacle. 

In Section ^ a theory of the surface waves generated by diffraction, that makes use of these 
geometrical properties, is presented. Throughout this work we keep, as a typical example, the 
diffraction of the light by a convex and opaque object, and, for the sake of simplicity, the light is 
represented as a scalar. By using the same method, analogous results can be obtained in sound 



diffraction (see the spectacular examples of creeping waves in acoustic diffraction in Ref. 14]) and, 
presumably, also in the diffraction of nuclear particles |15|| . 

A considerable part of the results that we obtain can be proved in a quite general geometrical 
setting: this is the case, for instance, of the proof of the existence of the diffracted rays, which 
derives solely from the non-uniqueness of the Cauchy problem. But, for other results, we have to 
restrict the class of obstacles to the surfaces of Besse type: i.e. the manifolds all of whose geodesies 
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are closed |Lq| . In particular, the sphere is a Besse surface. In this case all the results obtained 
by geometrical methods can be easily compared with those obtained by using standard methods 
based on the expansion of the amplitude in series (see subsection j3lj| ). When the obstacle radius 
is large, compared to the wavelength, the series converge very slowly, and the standard procedure 



suggests the use of the Watson resummation [17, 18]. However, these formal manipulations do not 



shed light on the actual physical process. Then, the geometrical approach is expected to be very 
useful in the investigation of more refined features, like ripples |19[| , which are very sensitive to 
initial conditions and size parameters. 

Finally, notice that our analysis will be limited to the geometrical theory of surface waves, 
whose effects are dominant in a small backward angular region, as it will be explained in subsection 



3.3 . The reader interested to a detailed analysis of which effect is dominant in the various angular 
domains, and, accordingly, to a systematic discussion of the transition regions, is referred to Ref. 
fi!| (in particular, Chapter 7 and Fig. 7.7). 

2 Riemannian Geometrical Optics 

2.1 Non-Uniqueness of Cauchy Problem in the Riemannian Manifold with 
Boundary: the Diffracted Rays 

In the variational derivation of geometrical optics, in particular for the laws of reflection and 
refraction, use is made of the Fermat's principle: the paths of the reflected or refracted rays are 
stationary in the class of all the paths that touch the boundary between two media at one point, 
assumed to be an interior boundary point. To introduce the paths of the diffracted rays it is required 
a generalization of the Fermat's principle, extended to include points, as well as arcs, lying on the 
boundary [19]. In our analysis, instead of the Fermat's, we use the Jacobi form of the principle 



of least action [20] which concerns with the path of the system point rather than with its motion 
in time. More precisely, the Jacobi principle states: if there are no forces acting on the body, 
then the system point travels along the shortest path length in the configuration space. Moving 
from mechanics to optics, Riemannian geometrical optics can be rested on the Jacobi principle, 
formulated as follows: the light rays travel along geodesies. 

In this context the diffraction by a convex, smooth and opaque object can be reconsidered as 
a Riemannian obstacle problem: the object is regarded as an obstacle which a geodesic can bend 
around, or which a geodesic can end at. Let IK denote the obstacle which is embedded in a complete 
n-dimensional Riemannian manifold (H, g), where g is the metric of HI and n > 3. Let us introduce 
the space M = H\K, where the obstacle IK is an open connected subset of H, with regular boundary 
dK and compact closure IK = IK U <9IK. Although most of the results illustrated below hold true in 
a very general setting, in the following we keep very often, as a typical example, HI = R 3 endowed 
with the euclidean metric. Finally, we are led to consider the space M* = M U <9IK (= H\IK), that 
has the structure of a manifold with boundary. 
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Now, two kinds of difficulties arise: the first one concerns geodesic completeness, that is the 
possibility to extend every geodesic infinitely and in a unique way; this uniqueness is indeed missing 
in M* at the points of the boundary. The second difficulty is related to the necessity of finding 
suitable coordinates at the boundary which allow for using the ordinary tools of the differential 
geometry, e.g. for writing the equations of the geodesies. 

Concerning the first point, the lack of geodesic completeness can be treated by introducing the 
notion of geodesic terminal (see Ref. [^TJ) to represent a point where a geodesic stops. Following 



Plaut |21 ], it can be proved that M* is the completion of M by observing that the set I of the 
geodesic terminals is nowhere dense and dM* = <9K = X. 

Regarding the second point, it is convenient to model the manifold with boundary on the half- 
space = {(a?i, x n ) G M n |x n > 0} C R n , where Mq" 1 denotes the boundary x n = of 
Thus, for a manifold with boundary, there exists an atlas {U a }, with local coordinates (x^ , ...,x"), 
such that in any chart we have the strict inequality x" > at the interior points, and = at the 
boundary points. The set <9M* of the boundary points is a smooth manifold of dimension (n— 1). 

In a Riemannian manifold without boundary, the geodesies are the solutions of a system of 
differential equations with suitable Cauchy conditions. The main result is the theorem which 
guarantees existence, uniqueness and smoothness of the solution for the Cauchy problem of this 
system. The case of Riemannian manifold with boundary is different from the classical one for 
what concerns uniqueness and smoothness of the solutions. 

In order to write the equation of a geodesic 7 of M*, we introduce suitable coordinates Xi 
(i = 1, . . . , re) adapted to the boundary, called geodesic boundary coordinates H ||, with x n defined 
as the distance from the boundary <9M*; then starting with arbitrary coordinates Xj, i < n on 
<9M*, we extend them to be constant on ordinary geodesies normal to cM* . Let denote the 
Christoffel symbols of the Levi Civita connection of M*, and x be the normal curvature of cM* in 
the direction of 7. Then, we have: 

X = ^ ~] %i %j ^ijn- (1) 

i,j<n 

When 7 is not in dM* the same expression occurs in the differential equations for 7, so that if we 
define x = off the boundary segments, then the differential equations can be modified to cover, 
in an integral sense, all points of 7, as follows [||, [|: 

&k = — ^ Xj Xj r^-fc, (k < n), (2) 



X ^ 1 %i x j r^jn- (3) 

i,j<n 

Equations (§, ||) can be implemented with the Cauchy conditions by providing the values of 
x(i) and x(t) (x = (xi, x n )) when t = 0: 

Xi (0) = b u (i<n), x n (0) = 0, (4) 
Xi(0) = Vi, (i < n), x n (0) = 0, (5) 
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where (61, ■■■,b n ) are the coordinates of the point b belonging to the boundary, and (v±, ...,v n ) are 
the components of the tangent vector at b. 
Now, let us consider the following cases: 

a) x n = 0. In view of equalities (|l|) and ([|), we remain with system (0) composed by (n — 1) 
differential equations that describe a geodesic on the (n — 1) dimensional smooth manifold 
<9M* (i.e. x n = 0), which is a manifold without boundary. In this case the standard theorem 
on the existence, uniqueness and smoothness of the geodesies satisfying the given Cauchy 
conditions holds true. We have a C°°-class geodesic which is completely contained on the 
boundary. 

b) x n (t) = if t < 0, and x n {t) > if t > 0. In view of the fact that x = off the boundary 
segments, eqs. (^, ||) describe, for t > 0, a geodesic segment belonging to the interior of 

o 

M* (i.e. M*). Moreover, eqs. (Q, |3|) describe, for t < 0, a geodesic segment belonging to 
the boundary. Notice that the r.h.s. of (g) is continuous, so that @ holds everywhere as 
it is. If x 7^ ° n boundary segments, from Eq. (||) it can be easily seen that x n fails to 
be continuous at the point b (t = 0). Nevertheless, in view of the Cauchy conditions (i.e. 
^n(0) = x n (0) = 0), the two geodesic segments (one belonging to the boundary, the other 
belonging to the interior of M*) glue necessarily at 6, and give rise to a geodesic of M* of 
class C 1 at the transitions point b. 

c) x n (t) = 0, if t > 0, x n (t) > if t < 0. The analysis is analogous to the previous one at point 
(b). We can say that two geodesic segments (one belonging to the interior of M*, the other 
one to the boundary) glue together at b and give rise to a geodesic of M* of class C 1 at the 
transition point b. 

d) x n (t) = if and only if t = 0. Following considerations strictly analogous to those developed 
above, we can say that eqs. (|2|, ^) describes two geodesic segments which belong to the interior 
of M*, and have only one point of contact with the boundary at b. Therefore, gluing at b the 
two geodesic segments, we obtain a geodesic of M* of class C°° . 

Since we supposed the obstacle to be compact with smooth and convex boundary, the analysis 
performed above is exhaustive. 

Summarizing, if x 7^ on the boundary segments, we can glue together, at the point b G 9M*, 
a geodesic segment belonging to the boundary with a geodesic segment belonging to the interior of 
M* (cases (b) and (c)), and, in these case, we have a geodesic of M*, which is of class C 1 at the 
transition point b. Thus, we have obtained the diffracted rays, which are precisely the geodesies 
of M* of this type. The situation considered in case (d) is of minor interest for our analysis, since 
those geodesies correspond to rays which are not diffracted even if they touch the obstacle. 

As a consequence of the non-uniqueness of the Cauchy problem at the boundary of the manifold, 
that can be made transparent by the use of equations (|2|) and (||), we can conclude by saying that 
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at each point of the boundary we have a bifurcation: the ray splits in two, one part continues as 
an ordinary ray (diffracted ray of class C 1 at the splitting point), the other part travels along the 
surface as a geodesic of the boundary. 

Finally, let us observe that the contact point b on the boundary is necessarily an elliptic point 
in view of the assumptions of convexity and compactness made about the obstacle. 



2.2 Conjugate Points (Axial Caustic) and Morse theorem in a Riemannian 
Manifold with Boundary 

First, let us consider a complete Riemannian manifold M without boundary. Given a geodesic 
7 = 7(i), < t < 1, consider a geodesic variation of 7, that is a one-parameter family of geodesies 
7s = 7s(t), (— e < s < e) such that 70 = 7. For each fixed s, 7 s (t) describes a geodesic when t varies 
from to 1. Each variation gives rise to an infinitesimal variation, that is, a certain vector field 
defined along 7. The points 7(0) and 7(1) of 7 are said to be conjugate if there is a variation 7 S 
which induces an infinitesimal variation vanishing at t = and t = 1 (see Ref. J22j). 

Set f(s,t) = "f s (t), and denote by D/dt and D/ds the covariant differentiation with respect to 
t and s, respectively. Denote by R the curvature transformation determined by df /dt and df/ds, 

\dt , 'bl) ~&t ~ \dldi~ dids~) ~dl' ( ' 

A vector field Y is a Jacobi field along 7 if satisfies the following differential equation: 

j Y = Y " + R( 1 ',Y) 1 ' =0, (7) 

where 7' = dj/dt, Y" is the second covariant derivative of Y along 7, and R (7', Y) is the curvature 
transformation defined by equality @. Two points p, q G M are conjugate along 7 if there exists a 
non-trivial Jacobi field Y along 7 such that Y(p) = Y(q) = 0. Finally, the multiplicity of the pair of 
conjugate points p, q £ 7 is given by the dimension of the vector space of the linearly independent 
Jacobi fields, along 7, that vanish at p and q. 

Let us now introduce the infinite dimensional space £l(M;p, q) of piecewise differentiable paths 
c connecting on M the point p with the point q: i.e. let p and q be two fixed points on M, and 
c : [0, 1] -» M be a piecewise differentiable path such that c (0) = p, c (1) = q. To any element 
c S 0(M;p, q) we associate an infinite dimensional vector space T c f2, which can be identified with 
the space tangent to Q in a point c. More precisely T c f2 is the vector space composed by all those 
fields of piecewise differentiable vectors V along the path c, such that V(0) = 0, ^(1) = 0. Next, 
we introduce the energy-functional 

2 



E(c) 



dc 
dt 



dt = gij x l x 3 dt, (8) 



where is the metric of M, and (x , x n ) are local coordinates on M. From the first variation of 
functional E we can deduce that the extremals of functional E{c) are represented by the geodesies 
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c(t) = 7(i), parametrized by t. Let us now consider the second variation E^ of the functional E 
along the geodesic 7 (it will be called the hessian of E), and let A be the index of the hessian, that 
is the largest dimension of the subspace of T 7 fi on which E^ is definite negative. One can then 
state the following Theorem: 

Theorem 1 (Morse Index Theorem Q) The index A of the hessian E** equals the number of 
points belonging to j(t) which are conjugate to the initial point p = 7(0), each one counted with its 
multiplicity. 

Now let us go back to our main concern, the Riemannian manifold with boundary M* . Following 
Alexander |23|] , we define as a Jacobi field any vector field Y along the geodesic 7 (whose boundary 
contact intervals have positive measure), which satisfies the following conditions: Y is continuous, 

o 

it is a Jacobi field of both M* and dM* along every segment of 7, and at each endpoint of a 
contact interval it satisfies the following equation: 

F (l r ) (tr) = p (i r ) (t ' +> ' (9) 

where P denotes orthogonal projections onto the hyperplane tangent to dM.* . In the Riemannian 
manifold with boundary we are forced to consider geodesies which lie both on and off dM*, and that 
are merely C 1 at the transition points (see the previous subsection). However, it is still possible to 
define a hessian form which is given simply by the sum of the classical formulae for dM* on contact 

o 

intervals, and for the interior of M* (i.e. M*) on interior intervals. Furthermore, we call 7(0) and 
7(1) conjugate if there is a non trivial Jacobi field Y along 7 whose limits at the endpoints vanish. 
Finally, the Morse index theorem can be extended to the Riemannian manifold with boundary. To 
this purpose, it is convenient to introduce the so-called regular geodesies in the sense of Alexander. 



Following Ref. 23], a geodesic 7 is regular if: 

a) all boundary-contact intervals of 7(f) have positive measure; 

b) the points of arrival of j(t) at dM* are not conjugate to the initial point 7(0). 
Then we can formulate the extended Morse index theorem as follows: 

Theorem 2 (Morse Index Theorem for Riemannian Manifolds with Boundary p3| ) Let 7 be a reg- 
ular geodesic; then the index A of E^ is finite and equal to the number of points belonging to j(t) 
conjugate to the initial point 7(0) (0 < t < 1), counted with their multiplicities. 

It still remains to check if the diffracted rays are regular geodesies. First of all we suppose that 

o 

the source of light is located in a point po SM* (point source); then the diffracted rays which we 

o 

are considering connect the point po and a point q SM*, placed at the exterior to the obstacle. 
If we assume that EI = R 3 (equipped with euclidean metric), then the diffracted rays consist of a 
straight line segment from po to the body surface, of a geodesic along the surface, and a straight 
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line segment from the body to q. Now, in view of the fact that these rays undergo diffraction, they 
cannot be simply tangent to the obstacle, but the contact-interval must be a geodesic segment of 
positive measure. Therefore condition (a) is satisfied. Concerning condition (b), we simply note 
that it is certainly satisfied in view of the fact that we consider obstacles formed by a convex and 
compact body embedded in M 3 , and the latter does not have conjugate points. 

2.3 Morse Index and Homotopy Classes of Diffracted Rays 

Hereafter the analysis will be limited to obstacles represented by manifolds all of whose geodesies 
are closed (i.e. Besse manifolds |l6|1 ). Merely for the sake of simplicity, henceforth we shall consider 
as obstacles spherical balls embedded inM 3 : i.e dK. = S 2 , equipped with standard metric. However, 
in view of the homotopic invariance, the main result of this section (see next ( ftO[ |TT|) ) holds true 
for any convex and compact Besse type manifold. 

o 

Consider the axis of symmetry A of the obstacle IK passing through the point po GM*, where it is 
located the source of light, and denote by A- the portion of this axis lying in the illuminated region 
(i.e. the same side of the light source), and by A + the portion of the axis lying in the shadow. The 
axis A equals A^ L)DL)A + , where D is the diameter of the sphere. Let us note that all the points of 
M* that do not belong to A + are connected to po by only one geodesic of minimal length, whereas 
the points q G A + are connected to po by a continuum of geodesies of minimal length that can 
be obtained as the intersections of M* with planes passing through the axis A. By rotating these 
planes, and keeping po and q G A + fixed, we obtain a variation vector field which is a Jacobi field 
vanishing at po and q. Thus, we can conclude that po and q G A + are conjugate with multiplicity 
one because the possible rotations are only along one direction. Then, in view of the extended form 
of the Morse Index Theorem, we can state that the index A of E** jumps by one when the geodesic 
7(i), whose initial point is 7(0) = po, crosses A + . 

Now let us focus our attention on the points q G M*\(A + U A^) (i.e. which do not lie on the 
axis of the obstacle, see Fig. ||). Consider the space X poq = (E 2 ) g n M*, where (E 2 ) g is the plane 
uniquely determined by the axis of the obstacle through po and the point q. X poq (which will be 
denoted hereafter simply by X) is an arcwise connected space, whose boundary in (E 2 ) g is a circle 
S which is a deformation retract of X. In view of these facts the fundamental group tti(X;pq) 
does not depend on the base point po and is isomorphic to 7Ti(5 1 ;l), where 1, described in a 
convenient system by the coordinates (0,1) (see Fig. ||), represents the contact point of A + with 
S 2 : tti(X;pq) ~ tti(S 1 ; 1) ~ Z. Let us consider, within the fundamental groupoid H\(X) of X, the 
set Hi(X;po,q) of paths in X connecting po with q, modulo homotopy with fixed end-points. We 
can then formulate the following statements: 

Proposition 1 Each element ofiri(X;po) is a homotopy class [a], with fixed endpoints, of a certain 
loop a : [0, 1] — ► X, in the space X , starting and ending at the point pq. 
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Figure 2: Geometry of the problem in the case of spherical obstacle. 70 is the unique geodesic 
of minimal length when the point q does not belong to the axial caustic (the z-axis). «o is a 
counterclockwise loop at po. 

Proposition 2 Each path cq from po to q £ M*\(A+. U A-), determines a one-to-one correspon- 
dence W between ir\{X]po) and the set T\-\{X]pQ,q). Such a correspondence can be constructed as: 
V [c] G Hi(X;po,q) : [c] — ► [c^Cq 1 ] G tti(X;pq), where the symbol denotes the concatenation of 
paths, and Cq 1 denotes the reverse path of cq. 

Proposition 3 i) In each homotopy class of iri(X;po) there is precisely one element of minimal 
length, ii) In each homotopy class of IL%(X;pQ,q), q G M*\( J 4 + Ui_), there is precisely one 
geodesic. 

Because of Propositions ^ and |3| we can establish a bijective correspondence w between geodesies 
from po to q (fixed) and integral numbers. Let us make precise our choices: 

a) q is fixed in X\(A+ U A_); we introduce in (E 2 ) g a reference system so that po belongs to 
the negative part of the z-axis, and the value of the coordinate x of the point q (i.e. x q ) is 
positive (see Fig. [2]); 

b) 70 is the geodesic from po to q of minimal length; 
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c) ao is a loop in X at po, such that: 

c') [ao] is a generator of iri(X;p ) (establishing an isomorphism tti(X;pq) ~ Z); 
c") ao turns in counterclockwise sense around the obstacle (see Fig. ||). 

Now, the correspondence w maps the geodesic 7 into the number 10(7) such that [7 * 7o~ 1 ] = 
10(7) [ao]. Note that w(-y) is the winding number of the loop [7 * Jq 1 ] determined by 7 (with 
respect to the chosen generator [ao]). 

We can also characterize the geodesies (from p$ to q) by a natural number 77.(7), which we call 
the crossing-number of 7, having a clearer geometric interpretation: 77(7) counts the number of 
crossings of 7 across the z-axis. For instance, 77(70) = 0. It is easy to see that the winding number 
determines the crossing- number through the following bijective correspondence: 

Z — ► N 

777 ► (2777 - 1) (777 > 0) (10) 

777 > -2777 (777 < 0) (11) 

so that our previous statement on the bijective correspondence between the geodesies and their 
crossing- number is established. 



3 Surface Waves in Diffractive Scattering 

3.1 Creeping Waves on the Sphere 

In this subsection we consider a Riemannian manifold without boundary MP; is a chart of MP, 
the local coordinates in $7 and (x) is the metric tensor. As usual g = |det (g%j)\ 
and = g~^ . Let us consider the Helmholtz equation 

A 2 u + k 2 u = 0, (12) 

where A2 is the Laplace-Beltrami operator which, for a function u G C°°(M n ), reads 



A 2 « = — — [>_.;" v 5— |« -LI) 

Now, we look for a solution of Eq. (|l^) of the following form: 

u(x,fc) = f A(x,(3) e ik ^ p) dp. (14) 



The principal contribution to u(x, k), as k — > +00, corresponds to the stationary points of in the 
neighbourhoods of which the exponential e lfc * ceases to oscillate rapidly. These stationary points 
can be obtained from the equation 5$(x, (3)/df3 = (provided that 9 2 $/5/3 2 7^ 0). Let us now 



E. De Micheli, G. Monti Bragadin, G.A. Viano 



12 



suppose that for each point x = (xi, ■~,x n ) <!> has only one stationary point; then the following 
asymptotic expansion of u, as k — > do, is valid |J: 

oo . 
m=0 V ' 

where /?o is the unique stationary point of $. The leading term of expansion (|l^) is: 

u(x,fc) =^ (x)e ifc * (x ' /3o) , (16) 

where 



d 2 $ 



dp 



7T 




exp Yi sgn {w)c 1 (17) 



//3=/3o 

By substituting the leading term (which will be written hereafter as Aexp{ik&}, omitting the 



index zero) into equation (12), collecting the powers of (ik), and, finally, equating to zero their 



coefficients, two equations are obtained: the eikonal (or Hamilton-Jacobi) equation 



i 9$ 9$ 
^ dxi dxj ^' 



and the transport equation 



whose physical meaning is the conservation of the current density. 

We can switch from wave to ray description by writing the eikonal equation in the form of a 
Hamilton-like system of differential equations by setting 

<93> 

" = K? (20) 
F = Ua'h.P-i), (21) 



and obtaining the characteristic system 

dxi 
Ik 



9 ij Pj, (22) 



dr dxi 2 ^ l ^ dxi ' 

where r is a running parameter along the ray emerging from the surface $ = const. 

Whenever the amplitude A becomes infinite, approximation (^) fails, and, consequently, such 
an approximation holds true only locally. Thus, we are faced with the problem of passing from a 
local to a global approximation, i.e. a solution in the whole space. The strategy for constructing 
a solution in the whole space consists in patching up local solutions by means of the so called 
Maslov-indexes in a way that will be illustrated below in the specific case of M n = S 2 . 
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Let the unit sphere be described by the angular coordinates 9 and </>. The matrix elements gij 
have the following values: gu = 1, 522 = sin 2 gyi = 521 = 0. In addition, we assume that the 
phase <3? and the amplitude A do not depend on the angle <§>. Then, eqs. (18) and (^) become, 
respectively 

1 f d f t _.___ nl ,2 d<S>\ 



| sin fl|^ =0 , e^rnr. (25) 



\smO\ [d9 \ 1 d9 J 

From Eq. ( g4j) we have: $ = ±9 + const. , and from Eq. (|25| ) we obtain: A(#) = const. 
I sin#| -1 / 2 , e 7^ n7r. Therefore, approximation (|l6|) becomes 



n(0,fc) = ^=e ± ^, M^, (26) 

V I sm0| 

where the terms exp{±i/c#} represent waves traveling counterclockwise (exp{ike}) or clockwise 
(exp{— ike}) around the unit sphere. Since the approximation (^) fails at e = nir, (n = 0, 1,2...), 
we have to consider the problem of patching these approximations when the surface rays cross the 
antipodal points = 0, tt, which are conjugate points in the sense of the Morse theory (see Section 



|2.3| ) . This difficulty will be overcome by the use of the Maslov construction [g4[ 25 , [10| . 

In order to illustrate the Maslov scheme in the case M n = S 2 , we reconsider the problem of the 
waves creeping around the unit sphere in a more general setting. The Hamilton-Jacobi equation 
(|l8l) is rewritten in the form: 

H = g i ip l p j =p 2 e + -± T -pl = l. (27) 
sin e 

In order to find the first integrals h of the system, we equal to zero the Poisson brackets: {h, H} = 0, 
to get: 

P<t> = ci, (28) 



,,2 
_]_ 

sin 2 



vl + -rh = 1 - ( 29 ) 



*±(M)= (P6d<t> + p d0) = c l <j>±$ 1 (9), (30) 



Furthermore, we have: 



(notice that pg = ±{1 — (c 2 / sin 2 0)} 1 / 2 ). Equation ( |29| ) defines a smooth curve in the domain 
< e < it, pg E M, which is diffeomorphic to the circle. The points where the tangent to the 
curve is vertical have coordinates (#o>0), (tt — 9q,0), with 9q = sin _1 (ci), c\ > 0. The cycles of 
singularities are given by the equations (see Ref. ||, Fig. 16): 

P4, = ci, p e = 0; 9 = e , 0<^<2vr, (31) 
P4> = ci, p e = 0; e = vr - 6 Q , < < 2vr. (32) 
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Then, the neighbourhoods of the points 9 = 9q and 9 = ir — 9q are badly projected on the con- 
figuration space (9,4>). However, the Maslov theory guarantees that it is possible to choose other 
local coordinates such that the mapping from the Lagrangian manifold to the selected coordinates 
is locally a diffeomorphism. In the present case it is easy to see that the mapping in the plane 
(pg,4>) is diffeomorphic (see Fig. 16 of Ref. ||). Once the local asymptotic solution in terms of 
(pg, 4>) has been computed, it is possible to return to the configuration space (9, </>) by transforming 
pg — > 9 through an inverse Fourier transform F^p -*0' ^ * s m deed the asymptotic evaluation (for 
large values of k) of F^* obtained again by the stationary point method, that gives rise to an 
additional phase-shift in the solution. In fact, the term exp{i(7r/4) sgn (<9 2 <3?/<9/3 2 )} of formula ( |i~7| ) 
is modified as follows. Instead of <£(#,</>), it must be considered the phase &(pg,4>) in the mixed 
representation coordinate-momentum; to have the phase described in the configuration space we 
move back from pg to 9 by means of the inverse Fourier transformation. Then, we have an expo- 
nential term of the form exp{ik[& + 9pg]} into an integral where 9 is regarded as a fixed parameter, 
whereas pg is the integration variable. Let $>' = <!> + 9pg. The stationary point is then determined 
by the equality 

8& <9<l 

^- = ^- + = 0, (33) 
dpg dpg 

and, therefore, it is given by that value of pg such that d&/dpg = —9. Next, we get: 

d 2 & d 2 4> 89 

op% dp l Q dpg 

Therefore, the additional phase factor exp{— i(ir/4) sgn {89 /dpg)} emerges. 

The negative inertial index of a symmetric non-degenerate (n x n) matrix A, called Inerdex, is 
the number of negative eigenvalues of the matrix || . The following relationship holds j|] : sgn A + 
2 Inerdex A = n (where sgn A is for the signature of the quadratic form associated with the matrix 
.A). In our case we have 

89 , 89 

— sgn — — = 2 Inerdex — 1, (35) 

ope 8pg 

[ 7T 89 1 ( 7F 89 it ~\ 

exp < — i — sgn — — > = exp < % — Inerdex — i— > . (36) 

I 4 dpg} [2 dpg 4 J 

We now focus our attention on the phase factor exp{i<5} = exp{i ^ Inerdex (J^)}, which is the 
relevant term in the analysis of the crossing through the critical points. In close analogy with the 
Morse index which gives the dimension of the subspace on which is negative definite, and whose 
jumps count the increase of this dimension (see Section [^), here it must be evaluated the number 
of transitions from positive to negative values of (89 /dpg), or vice- versa. With this in mind, we 
cover the curve: p 2 , + (c 2 / sin 2 9) = 1, diffeomorphic to a circle, with four charts: let U 1 be the chart 
which lies in the half-plane pg < 0, U s be the chart which lies in the half-plane pg > 0, and U 2 and 
U A be the charts which lie in the neighbourhoods of the points (#o>0) and (tt — 9q,0), respectively. 



then: 
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Next, the evaluation of (d9/dpo) gives 

96 cip e 



(0 < a < 1). (37) 



dpe (i-^/HW) 

Let us consider a path I, counterclockwise oriented, which connects two points r' and r" belonging 
to U 1 and U 3 , respectively. The Maslov index Ind (I), defined as the variation of Inerdex (dO/dpo) 
along I, reads 

/ QQ \ / QQ \ 

Ind(Z) = Inerdex I — — I — Inerdex I — — ] = —1, (38) 



\dpej r » \dpe 
since Inerdex {d6/dpe) r „ = 0, because pe > 0, while Inerdex (d6/dpo) r , = 1, because pe < 0. 
Accordingly, we have a phase- factor: exp{i<5} = exp{— in/ 2}. Conversely, for any other path /' 
connecting two points lying in the same chart U 3 (or IA ) that does not intersect the cycle of the 
singularities, we have Ind (/') = 0. 

Remark 1 Unfortunately there is a very unpleasant ambiguity concerning the sign of Ind (/) (see 
Ref. ^4|). Here we follow the Maslov prescription [25]: the transition through a critical (focal) 
point in the direction of decreasing (80/dpg) increases the index Ind(Z) by one; the transition in 
the direction of increasing (d9/dpg) decreases the index Ind(Z) by one. 

Now, let us consider the patchwork across the critical points for the derivation of the connection- 
formulae. As we have seen previously, u(k, 6, 4>) can be approximated, for large values of k, by the 
expression: Aex.p{ikQ±} (see formula (p0|)). After crossing a critical point, say {9q, 0), an additional 
— 7r/2 phase-shift arises, and, accordingly, the solution becomes Aexp{i[k&± — vr/2]}. By taking 
the limit c\ — ► 0, we obtain the following connection-formulae: 

a) For the rays traveling around the sphere in counterclockwise sense, we have: 

after the 
crossing 

e lke . e * P*-*], (39) 

(we use the convention of taking positive the arcs counterclockwise oriented). 

b) For the rays traveling around the sphere in clockwise sense, we have: 

after the 
crossing 

e~ ike . e-^-il, (40) 

(the arcs oriented in clockwise sense are taken negative). 

Finally, in view of the homotopic invariance of the Maslov index ||, the bijective correspondence 
between the winding number, associated with the fundamental group tti(X;pq), and the crossing- 
number (discussed in Section ^j^) can now be extended to the Maslov phase-shift as follows: 

Z — ► N — * Maslov phase-shift 

m — ► (2m - 1) — > - | (2m - 1) (m > 0) (41) 
m — ► - 2m — ► | (-2m) (m < 0) (42) 
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In particular, notice that, for every complete tour around the sphere, both the counterclockwise 
and the clockwise oriented rays acquire a factor (—1), due to the product of two Maslov factors. 

Remark 2 i) In the literature the term creeping-waves usually indicates waves creeping along the 
boundary and continuously sheding energy into the surrounding space. This is also the meaning of 
this term in the present paper. However, only in the present subsection, this term has been used 
with a slight different meaning: the waves are creeping around the obstacle, but they are supposed 
not to irradiate around. In fact, in this subsection we have considered a manifold without boundary. 
In the next subsection the problem will be reconsidered in its full generality as an obstacle problem 
in a Riemannian manifold with boundary, and we shall evaluate the damping factors associated 
with the rays that propagate along the boundary, and that shed energy into the surrounding space. 



ii) Guillemin and Sternberg, in their excellent book [26] on Geometric Asymptotics, give an analysis 
of Maslov's indexes, and illustrate the related application in a very general setting, with particular 
attention to the geometric quantization. They also calculate the ir/2 phase-shift at the crossing of 
the caustic by the use of the Morse theory. 

3.2 Airy Approximation and Damping Factors 

The diffraction problem concerns with the determination of a solution u s (x, k) of the reduced 
wave equation (|l^) in the exterior of a closed surface S that, in our case, is a sphere of radius R 
(embedded in R 3 ). An incident field Uj(x, k) that satisfies Eq. (^) is prescribed, and « s (x, k) is 
also required to satisfy the radiation condition 

lim |x| ( -f-rM s - iku s J = 0, (43) 

|x|^oo V^l x l / 

(where |x| = Q^iLi x 1 ) 1//2 j x i being the cartesian coordinates of the ambient space M 3 ). Finally, 
the total field u(x, k) must also satisfy a boundary condition on the obstacle (see below formula 

In order to solve this problem we intend to apply once again the method of the stationary phase 



to an integral of the form (14). But, in this case, the surface S (boundary of a Riemannian manifold 
M*) is the envelope of the diffracted rays: it is a caustic. We shall prove in the next subsection that 
the approximation given by eqs. ( |18| ) and ( |19D fails locally on this domain, since the amplitude 
becomes infinite. Each point P outside the caustic lies on the intersection of two diffracted rays 
which are tangent to the border of the diffracting ball, and crossing orthogonally two surfaces of 
constant phase ^ (for simplicity, we are now considering rays that are not bending around the 
obstacle). When the point P is pushed on the boundary, the curves of constant phase meet forming 
a cusp, and, accordingly, the stationary points f3\ and fa of the two phase functions ^ coalesce. 



In such a situation the standard method, used in Section 3.1 fails, and a different strategy must be 



looked for. An appropriate procedure is the one suggested by Chester, Friedmann and Ursell [11], 
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that consists in bringing the phase function $ into a more convenient form by a suitable change of 
the integration variable £ <-> /3, implicitly defined by 

$(x,/3) = e (x) + F(x,0. (44) 

After this change, integral (14) reads 

u(x,k) = e ifce ° (x) J s(x,0 e ikF{x >V (45) 

This expression is indeed similar to the integral (|i~4|), with the phase- function replaced by F, and 
with an additional oscillatory factor exp{i/c0o} in front. 

Now, we have two stationary points /?i(x) and /?2( x ) that coalesce, and our goal is to choose a 
transformation such that to these points there correspond the points where dF/d£ vanishes. This 
result can be achieved by setting 

F(x,0=-^ 3 + Po(x)e (46) 



In fact, dF/d^ = — £ 2 + /?o( x ) = gives £ = ±^fp Q . Then, from equalities (44) and ( (46| ) we obtain 
the following relationships: 

x) = |pg /2 (x) + G (x), (47) 
^(/3 2 ,x) = ~p 3 /2 (x) + 8 (x), (48) 



that yield 



9 (x) = i {$(/?!, x) + cD(/3 2 ,x)}, (49) 
|p5 /2 (x) = i{$(A,x)-*(A,x)}. (50) 



In the case /3j = we have /3o(x) = and ©o(x) = $(/3i,x) = <I>(/?2, X ). If equalities (49, |5C 
are satisfied, the transformation £ <-> (3 is uniformly regular and 1 — 1 near £ = (see Ref. PI). 
From eqs. ([45|) and fl46| ) it follows that the most significant terms in the expression of it(x, &), for 
large fc, can be written in terms of the Airy function Ai(-) and of its derivative Ai'(-). Such a 
representation of u(x, k) led Ludwig |0| to propose the following Ansatz: 

u(x,fc) = e ite »W { 5o (x)Ai(-A: 2 / 3 po) + ^- 1 / 3 / lo (x)Ai / (-A: 2 / 3 po)} , (51) 

where <7o( x ) an d ^o( x ) are respectively the first terms of the two formal asymptotic series: 
Z™=o9m(x)k- m , and Z™= hm(x)k- m . 

Unfortunately, representation ( |5l|) cannot be made to satisfy the boundary condition at the 
surface of the obstacle. In fact, on any point of the surface S the two stationary points coalesce, 
and the phase of the diffracted ray can be assumed to coincide with the phase of the incident ray, 
that is x) = $(/?2,x) = 3>j(x) = phase of the incident wave. Consequently, po( x ) is identically 
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zero on S (see Eq. (|50|)). It is therefore necessary to introduce an appropriate modification of the 
asymptotic representation fl5l|) in order to satisfy a boundary condition which, in its general form, 



reads p7| : 

^- + ik 2 / 3 (u = (onS), (52) 

where £ is a smooth impedence function defined on S, and du/dN is the normal derivative. For 
C = 0, and C = oo the relationship ( |52"| ) reduces respectively to the Neumann's du/dN = 0, and to 
the Dirichlet's u = boundary condition on S. Following Lewis, Bleinstein and Ludwig p7p , we 
introduce the new Ansatz, which is a slight modification of the asymptotic representation (|5l|): 

u(x, jfe) = e ifc0(x) { 5 (x)A(-fe 2/3 p) + ifc- 1/3 /i(x)A'(-fc 2 / 3 p)} , (53) 

where: 

^( x ) ~ Z^-^TT' ~ 2^ -^7T' ( 54 ) 

m=0 m=0 

p(x) = /O0 (x) + A:- 2 / 3 /Ol (x), (po(x) = on S), (55) 
G(x) = 9o(x) + A:- 2 / 3 G 1 (x), (56) 
A{t) = Ai(te 42 "/ 3 ). (57) 

Now, by inserting formula ( |53| ) into the reduced wave equation (12), and collecting the coefficients 
of k mj,3 Aexp{ikQ} and k m / 3 A' exp{i/cB}, we obtain the following set of equations: 

(VG ) 2 + po(Vp ) 2 = l, (58) 

V6 • Vp = 0, (59) 

2V6 • V0i + Pi{V Po ) 2 + 2p V Po ■ Vpi = 0, (60) 

V6 • Vpi + V6i • Vp = 0. (61) 



Then the Ansatz (|5^ ) can be inserted into the boundary condition (|52|). Here we limit ourselves 
to report the result which is relevant in our case considering, for simplicity, the Dirichlet boundary 
condition, i.e. putting £ = oo in formula (|52|). The reader interested to the details of the calculations 
is referred to Ref. j27j. It can be easily seen from Eqs. ( |60|) and ( |6T|) that &\ acquires an imaginary 
term, that depends on the curvature of the obstacle. If the obstacle is a sphere of radius R, we 
have: 

... 1 / o \ 2 / 3 r l 

0W = const. - - qi e-/ 3 - dr, (62) 



2" \RJ 

I being the length of the arc on the surface of the sphere described by the diffracted rays, and % 
denoting the i-th zero of the Airy function. For each qi one gets a solution ( |53l) of the reduced wave 
equation with boundary condition (|52|), and, correspondingly, an infinite set of damping -factors aij 
given by 

on = e" ft/ , fa = const. |$| k l/3 R~ 2/3 , (63) 
which depend on the curvature of the obstacle. 



E. De Micheli, G. Monti Bragadin, G.A. Viano 



19 



Remark 3 In connection with the results of this subsection we want to mention the deep and 
extensive works of Berry & Howls [28, 29], and Connor and collaborators [ 30 1 on the asymptotic 
evaluation of oscillating integrals when two or more saddle points coalesce. Moreover, it is worth 



mentioning the paper by Berry [31|, where the first general statement of uniform approximations 
using catastrophe theory has been done (for the application to special cases see also Ref. ]32|]). 



x 




Figure 3: Diffr active scattering: geometry of the contribution of two grazing rays to the scattered 
amplitude. The ray below (solid line with open arrows) travels 9q radians in the counterclockwise 
direction along the boundary of the obstacle, and leaves the obstacle in the direction of the scattering 



angle 9 S (see formula (72)). The ray above (solid and dashed lines with filled arrows) travels in 
the clockwise direction along the surface of the obstacle and crosses the axial caustic twice before 
emerging in the direction of the scattering angle S (see formula (|73] ) ) . 



3.3 Surface Waves in Diffractive Scattering 

In the scattering theory the far field diffracted by the obstacle must be related to the incoming field 
whose source is located at great distance from the obstacle. Then, we have to consider the limit 
obtained when the source point po is pushed to — oo and the observer to +oo. Let us introduce 
an orthogonal system of axes (x, y, z) in R 3 whose origin coincides with the center of the sphere, 
and such that the z-axis, chosen parallel to the incoming beam of rays, is positively oriented in 
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the direction of the outgoing rays (see Fig. ||). Let us introduce, in addition, the coordinates on 
the sphere: ro is the radial vector, 0q is the azimuthal angle, Go is the angle measured along the 
meridian circle from the point of incidence of the ray on the sphere. Finally, let r be a parameter 
running along the ray; in particular we have (on the sphere): r = R6q (R being the radius of the 
sphere). Now, consider the rays that leave the surface of the sphere after diffraction. In view of the 
fact that the interior of the space, outside the obstacle, can be regarded as a Riemannian euclidean 
space without boundary, eqs. (22, 23) hold true, and, in this case, read 

I = p- m 
t = ° («> 

where r = [x, y, z), and p = V$ (3> = phase function). From eqs. (|o4|, |65| ) we get: 

r(6» o ,0o,r) = r ((p ,6 ) + rp o (0 o , Go), (66) 

where po is the unit vector tangent to the obstacle where the ray leaves the sphere. 
Let us now focus our attention on the ray hitting the sphere at the point of coordinates (-R, 0, 0), 
and then traveling in counterclockwise sense. The components of the radial vector are (see Fig. 
|^): ro = (— -Rcos Go cos <f>o, Rcos Go sin 4>o, -Rsin Go), and po = (sin6*o cos <j>o, — sin^o sin</>o, cos Go)- 
Substituting these expressions in Eq. (p6| ) we have: 

r(6 l o, 4>o,t) = (—Rcos Go cos 4>o + r sin#o cos (f>o,Rcos Go sin</>o — r sin Go sin 4>o, R sin Go + r cos #o)- 

(67) 

<9(x, y , z) 

Now, the domain where the Jacobian J = ' — ■ — - = t(R cos Go — Tsin^o) vanishes is composed 

5(6» ,(/>o,r) 

by: 

i) the surface of the sphere, where r = 0; 

ii) the semi- axis represented by r = f = R cot Gq- 



We can rewrite the transport equation (19) in the following form ||: 

i£<JA»)=0, (08) 

whose solution A = C/\f~J indicates once again that the amplitude becomes infinite for J = 0, i.e. 
on the caustic. In order to treat the scattering problem we must perform the limit for r — > +oo 



(r = |r|). Since r = \/r 2 — R 2 , then as r — > oo, r tends to r, and from the expression of the 
Jacobian, we obtain: y/j — > i r \/ sin #o . 

Now, let us revert to formula (^) in order to evaluate the contribution due to the diffracted 
ray. Recalling the asymptotic behaviour of Ai and Ai', and by noting that in the present geometry 
Go = -R#o + r ) an d that the length of the arc described by the diffracted ray on the sphere is 
/ = RGo, we obtain, for k — ► +oo and for large values of r (see formula (^|) and Ref. p7f| ): 

. . p ikr p -f}iRQo p ikRQo 

Ui (k,G ,r) — > C d \k) — e ——L=^, (69) 
fc^oo r zvsmt/n 

r — ».qo v u 



E. De Micheli, G. Monti Bragadin, G.A. Viano 



21 



(0 < #o < 7r), where is the exponent of the i-th damping factor Qj obtained in Section 3J2. For 
what concerns the diffraction coefficients cj- d \k), they have been derived by Lewis, Bleinstein and 
Ludwig [^7| by assuming that the amplitude of the diffracted wave is proportional to the amplitude 
of the incident wave. In the present geometry these coefficients read 

^-—■iSW- (70) 

(qi being the zeros of the Airy function), and coincide with the expression derived by Levy-Keller 
(see Table I of Ref. ||). Since the roots of the Airy function are infinite, we have, correspondingly, 
a countably infinite set of modes Ui . For the moment we focus our attention on the damping factor 
ao, whose exponent (5q is the smallest one, and, accordingly, on the mode uq that, for large values 
of k and r, reads 

u (k,9 ,r) — > -W^Xk) (71) 

fc-+go r y/sm9 

(0<6 < vr, A = R(k + i0o)). 

In order to evaluate the contribution to the scattering amplitude, the scattering angle 9 S must be 
related to the surface angle 9q. To distinguish between the contribution of the counterclockwise rays 
and that of the clockwise ones, we add the superscript '+' to all what refers to the counterclockwise 
rays, and the superscript '— ' to all the symbols referring to the clockwise oriented rays[]. With this 
convention, and observing that 6 S = 9q (see Fig. ||), the contribution to the scattering amplitude 
/i j of the counterclockwise grazing ray that has not completed one tour around the obstacle^ can 
be written as 

/ ( + 0) (M.) = Cf(*) , (0 < 9 S < vr). (72) 

Notice that in formula ([72]) the factor exp{— m/2] corresponds to the Maslov phase-shift due to 
the fact that the ray crosses the axial caustic once (see Fig. ||). 

Analogously, the contribution to the scattering amplitude fZ Q ^ of a diffracted ray which travels 
around the sphere, in clockwise sense, without completing one tour around the obstacle (see Fig. 
|3|) can be evaluated. To this purpose, it is convenient to consider the asymptotic limit of |J| 1//2 , 
which is given by r | sin 9q \ 1//2 . Since the scattering angle 9 S is related to the surface angle #q~ as 
follows: 9q = 2tt — 9 s , and by noting that this (clockwise oriented) grazing ray crosses the axial 
caustic (i.e. the z-axis) two times (see Fig. |3|), we have: 

fr Q M s ) = (-1) C^(k) , (0 < 9 S < vr), (73) 



smi 



1 In the present geometry the simple convention adopted in Section |3.l| is ambiguous and cannot be used. In the 
present case the rays turning in clockwise sense hit the sphere at a point which is antipodal with respect to the point 
where the rays, turning in counterclockwise sense, hit the obstacle. 

2 In the notation Q , , the first zero of the subscript indicates that the grazing ray has not completed one tour, 
whereas the second one indicates that we are taking the smallest exponent /3o- 
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where the factor (—1) is precisely given by the product of two Maslov factors^ (see Section |3.1|) . 
Adding /i ^ to fZ Q y we obtain the scattering amplitude /(o,o)(&> @s) due to the rays not completing 



one tour around the obstacle: 

f m (k, 0s) = / ( + 0) + / ( o, ) = -* Ci d \k) == , (0 < 9 S < tt). (74) 



ysini _ 

Now, we are ready to take into account the contribution of all those rays which are orbiting 
around the sphere several times. Let us consider the rays describing n (n £ N) tours around the 
obstacle. Since the surface angles 9^ n are related to the scattering angle 9 S as 9$ n = 9 S + 2im, 
9q h = 2tt -9 s + 2irn (n = 0,1,2,.. .), we have for < 9 S < tt: 



The factor (—1), at each n, is due to the product of two Maslov phase factors, corresponding to the 
fact that both the counterclockwise and the clockwise rays cross the z-axis (i.e. the axial caustic) 
twice for each tour (see Section [O]). 
Now, we exploit the following expansion: 

1 oo 

= e wAo yV^ngfitanAo (Im Xq > q\ ^ 

2 cos 7rAn ' 
n=0 



By the use of formula (76), the amplitude (|75| ) can be rewritten as 

/ (fc, S ) = -C#°(fc) e-/ 4 -Z , (0 < S < tt). (77) 

(z cos 7rAo) vsmc/ s 

The r.h.s. of formula (|77| ) contains the asymptotic behaviour, for |Aq| — > oo and | Aq | (tt — 9 S ) >> 1, 



of the Legendre function P Aq _i (— cos S ) times y / 27rAo |3^]. Then, writing -f\ _ i (— cos 9 S ) in place 
of its asymptotic behaviour, we have for |Ao| — » oo (0 < # s < 7r): 

/a v 7 27rAn -P\ i(— COS0 S ) 

Hk,9 s ) * -Ci d \k)e^ ° -, (78) 

2 COS 7TAo 

that, by putting /xq = Ao — 1/2, becomes 

/o(M.) * cf (&) e-/ 4 ^7(^Tl) ^ o( ~ COSgs) . (79) 

2 sin 7T/uo 

Now, if we consider the countable infinity of damping factors 04 (see Section j3^2j ), we obtain an 

infinite set of creeping waves, whose angular distribution is described by the Legendre functions 

P Mi (-cos0 s ), where m = Aj - 1/2, Aj = R(k + i/3i). 

Let us recall that the simplest approach to diffraction scattering by a sphere is the expansion 

of the scattering amplitude in terms of spherical functions. But, when the radius of the sphere 

3 It is easy to prove, through arguments based on symmetry, that the negative z-axis is a caustic for the rays 
propagating in the backward emisphere. 



E. De Micheli, G. Monti Bragadin, G.A. Viano 



23 



is large compared to the wavelength, these series converge so slowly that they become practically 
useless (see Refs. [17, f\M). A typical example of this difficulty is given by the diffraction of 



the radio waves around the earth. In order to remedy this drawback, Watson [18[] proposed a 
resummation of the series that makes use of the analytic continuation from integer values I of the 
angular momentum to complex A- values (or /z- values in our case). In this method the sum over 
integral I is substituted by a sum over an infinite set of poles corresponding to the infinite set of 
creeping waves, whose angular distribution is described by the Legendre functions Pm(— cos 8 S ). 
Let us note that at 9 S = 0, cos# s ) presents a logarithmic singularity [[l?]], and, consequently, 



approximation (7£) fails. Furthermore, at small angles the surface waves describe a very small 
arc of circumference, and the damping factors on are close to 1; therefore, we are obliged to take 
into account the contribution of the whole set of creeping waves. On the other hand, at 9 S = ir, 
Pa. (— cos 6 S ) = 1 and, furthermore, at large angles, the main contribution comes from that creeping 
wave whose damping factor ao = exp(— (3qI) has the smallest exponent /3q. Therefore, in the 
backward angular region the surface wave contribution is dominant, and the scattering amplitude 
can be approximately represented as follows: 

f(k,0 s )~fo(k,e s ) = g o (k)P flo (-co S 9 s ), (0 S >O), (80) 

where G (k) = C^ d) (k)e^ (2/x + I) 1 ' 2 / smvuo- 



4 Conclusions 

Let us conclude with a brief remark on the difference between complex and diffracted rays. Complex 
rays are well known in optics in total reflection, where they describe the exponentially damped 
penetration into the rarer medium associated with surface waves traveling along the boundary. 
Then, it is very tempting to describe diffraction in terms of complex rays. Conversely, in the 
theory of diffraction which we present in this paper, we do not make any use of complex rays, 
but rather we introduce the diffracted rays, that are real. The first difficulty which emerges in this 
approach is the proof of the existence of these diffracted rays. This problem is solved by showing the 
non-uniqueness of the Cauchy problem for the geodesies in a Riemannian manifold with boundary. 
Then, since the border of the obstacle is the envelope of the diffracted rays, and the standard 
method of the stationary phase cannot be used, we are forced to apply a modified stationary phase 
method due to Chester, Friedman and Ursell [11]. By using this method we derive an infinite set 



of damping factors associated with the waves creeping around the obstacle. These damping factors 
have again a geometrical nature, since they depend on the curvature of the obstacle. In conclusion, 
it worth remarking that if we pass from optics to mechanics, and we consider particle trajectories 
instead of light rays, the splitting of geodesies at the boundary introduces a probabilistic aspect 
reflecting the fact that, at any point of the boundary, the particle can continue to orbit around the 
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obstacle or can leave the obstacle itself. Therefore, using the damping factors is a way to connect 
probability, proper of semi-classical mechanics, to geometry: i.e. the curvature of the obstacle. 
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